Non-equilibrium entangled steady state of two independent two-level systems 
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We determine and study the steady state of two independent two-level systems weakly coupled to 
a stationary non-equilibrium environment. Whereas this bipartite state is necessarily uncorrelated 
if the splitting energies of the two-level systems are different from each other, it can be entangled if 
they are equal. For identical two-level systems interacting with two bosonic heat baths at different 
temperatures, we discuss the influence of the baths temperatures and coupling parameters on their 
entanglement. Geometric properties, such as the baths dimensionalities and the distance between 
the two-level systems, are relevant. A regime is found where the steady state is a statistical mixture 
of the product ground state and of the entangled singlet state with respective weights 2/3 and 1/3. 

PACS numbers: 03.67.Bg,03.65.Yz,05.70.Ln 



I. INTRODUCTION 



For a quantum system, the influence of the surround- 
ings plays a role at a fundamental level. When the envi- 
ronment is taken into consideration, the system dynam- 
ics can no longer be described in terms of pure quantum 
states and unitary evolution. An open quantum system 
is generally in a statistical mixture of pure states. This 
has an important consequence for multipartite systems. 
As is well known, correlations between quantum systems 
cannot be completely understood in classical terms 
There exist states which are not classically correlated 
and lead to correlations with no classical counterpart, 
as clearly shown by violations of Bell inequalities for in- 
stance 0- They are said to be entangled. Whereas al- 
most all pure states are entangled, this is not the case 
for mixed states. In the space of mixed states, the set 
of non-entangled, or separable, states has a finite volume 
Q. An interesting consequence of the geometrical prop- 
erties of this set is that the state of a multipartite open 
system can be entangled for finite periods of time, in the 
course of its evolution, and separable at infinite time or 
vice versa 

The most common environment is a heat reservoir. If 
the considered system is weakly coupled to an infinite 
number of degrees of freedom, initially in thermal equi- 
librium, it relaxes, in general, to a thermal state with 
the temperature of its surroundings. In such an environ- 
ment, it is clear that, in the absence of direct interac- 
tions between the subsystems of a multipartite system, 
these subsystems are uncorrelated at long times. In other 
words, any initial correlation, quantum or classical, be- 
tween independent subsystems is generically destroyed 
by a thermal bath. Moreover, for the geometric reasons 
mentioned above, quantum disentanglement can occur 
in a finite time [5|, |6|. Furthermore, the disentangling 
influence of the environment also exists when no energy 
is exchanged between the system and its surroundings, 
whereas, in this particular case, classical correlations can 
persist 



However, when independent systems interact with a 
common environment, the indirect interaction between 
them, mediated by this environment, may have a positive 
impact on their entanglement. Recent results evidence 
the existence of this influence. It has been shown that 
a transient entanglement, between initially uncorrelated 
systems, can be induced by a thermal bath, for both non- 
dissipative @ and dissipative [Iol - [l2j couplings. It has 
also been obtained that, in the limit of infinitely close 
non-interacting systems, some special entangled states 
are not affected by the environment [l3j . In this limiting 
case, the considered multipartite open system has not a 
unique steady state, which is exceptional, and hence the 
entanglement evolution depends on the system's initial 
state. 

In the above cited dynamical studies, the environment 
is in thermal equilibrium and thus a relaxation dynamics 
towards a unique steady state necessarily means decay of 
correlations, both quantum and classical, between non- 
interacting systems. This may not be the case for a non- 
equilibrium surroundings. Stationary entanglement has 
been found in the presence of particle [lil . [TEj or energy 
flow [H, [13 ■ However, in these studies, entanglement 
occurs between systems that interact with each other di- 
rectly, via a two-level system, or via strong coupling to a 
heat bath, and this interaction plays an essential part in 
the development of entanglement. Such a strong interac- 
tion has been shown to be unnecessary for a different kind 
of non-equilibrium environment [l8j . In the presence of a 
classical oscillating field, the steady state of two two-level 
atoms, interacting with each other only via weak coupling 
to the electromagnetic vacuum, can be entangled. 

In this paper, we consider two independent two-level 
systems (TLS) coupled to a steady non-equilibrium envi- 
ronment. Examples of such surroundings are illustrated 
in FigJT] They consist of several heat baths at different 
temperatures. These are not the only possible examples 
and the two following sections are relevant to other en- 
vironments. In section mi we present the model used to 
describe two non-interacting TLS in a stationary environ- 
ment. In section UTTl we first study the steady state of a 
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FIG. 1: Schematic representation of non-equilibrium environ- 
ments of two independent TLS. The depicted environments 
consist of several heat baths at different temperatures. 



general system weakly coupled to its surroundings, and 
we then apply our approach to the case of a system con- 
sisting of two independent TLS. The system steady state 
is obtained, in the weak coupling limit, by solving per- 
turbatively an eigenvalue problem, which is derived from 
the system dynamics for arbitrary coupling strength. As 
far as one is interested only in the stationary state, no 
other approximation, such as a Markovian assumption, 
or elaborate method, such as a projection superoperator 
technique, are needed [III H3| • In section IIV1 we focus 
on the special case of an environment that consists of 
bosonic heat baths at different temperatures. It is shown 
that two baths are enough to induce stationnary entan- 
glement of two identical TLS. The influence of the two 
baths temperatures and of the coupling parameters is dis- 
cussed in some detail. Finally, we summarize our results 
in the last section. 



to lead to a decohering and dissipative reduced dynamics 
of the TLS. 

As the initial state of the complete system, we consider 



Q = J2 r ki\k)(l\®P£ 



(3) 



k.l 



where pe commutes with He . The two-TLS system and £ 
are initially uncorrelated. As we will see below, the condi- 
tion [pe, He) = implies the stationarity of relevant cor- 
relation functions of £ . Typical environments we are in- 
terested in are made up of several heat baths at different 
temperatures T n , as sketched in FigJT] In this case, the 
environment Hamiltonian and initial state read, respec- 
tively, as Hg = ^2 Hg n where n runs over the heat reser- 
voirs and [H en , H £ „>] = 0, and p £ oc n„ exp(-H Sn /T n ), 
and commute with each other. Throughout this paper, 
we use units in which h = fc^ = I . 



III. NON-EQUILIBRIUM STEADY STATE 

In this section, we first derive a matrix equation for the 
steady state of a generic open system S initially uncor- 
related with its environment £. More explicit equations 
are then obtained for a steady environment and in the 
limit of weak coupling between S and £ . This weak cou- 
pling approach is applied to the two TLS described by 
the Hamiltonian (p} . In this case, the steady state equa- 
tion can be solved. The result is radically different for 
Ai 7^ A 2 and A : = A 2 . 



General case 



II. MODEL 

The total Hamiltonian of two independent TLS and 
their environment £ can be written as 



(i) 



t (*) 
w'a_ 



H £ (1) 



= v] 



where A, are the level splittings of the TLS, v. 
and Wi are operators of £ and He is the self-Hamiltonian 
of £. The Pauli operator a z has eigenvalues ±1 and 
the corresponding eigenstates are denoted by |±) 



(*) 

operators <r^ 



then read a 



The 
We 



introduce, for further use, the following notations : 

|I) = |+}l| + >2 , |2> = |+>l|->2, 

|3) = |-}x|+) 2 , |4) = |-> 1 |-> 2 . 



(2) 



Two TLS interacting with their environment but not 
directly with each other can always be described by a 
Hamiltonian of the form (|T|). The system £ is assumed 
to consist of an infinite number of degrees of freedom and 



In general, under the influence of its environment £ , 
a system S relaxes to a steady state determined by its 
self-Hamiltonian and by its interaction with £ . If £ is 
in thermal equilibrium and S interacts with it weakly, 
this state does not depend on any detail of the intrinsic 
dynamics of £ or of the coupling between S and £ . But, 
as we will see, this is a very particular case. To determine 
the steady state of S, we first write its reduced density 
matrix, at positive times t, as 



Z7T 



dze- lzt TT f 



R+ir) 



(4) 



where Trg denotes the partial trace over £ , r\ is a positive 
real number, and f2 is the initial state of the total system 
S + £ . The Liouvillian C is defined by C . . . = [H, . . .] 
where H is the Hamiltonian of S + £ . This Hamiltonian 
can be decomposed as H = H$ + Hi nt + He where Hg 
and He are the self-Hamiltonians of S and £, respec- 
tively, and Hi n t accounts for the interaction between S 
and £ . The condition TTe(peHi nt ) — can be assumed 
without loss of generality. It can always be satisfied by 
appropriately redefining H$ and H int . The eigenstates 
and eigenenergies of Hg will be denoted by \k) and in 
the following. 
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1. Steady state equation 

For an initial state ft of the form ([3]), the matrix ele- 
ments ff.i(z) = (k\Ti£[(z~ £) -1 f2]|/) of the Laplace trans- 
form of p, are given by 



I'M 



where the functions 



:l( z ) = E Fkl,k>l>(z) 



Tk'V 



(5) 



k'l' 



Tki,k>i>{z) = (fc|Tr £ [(z-C) 1 \k')(l'\®pe]\l) (6) 

depend only on the environment part of the initial state 
(U). Equation (JS|) can be read as a matrix relation be- 
tween two column vectors r and f (z) with elements r k i 
and rki(z), respectively, and a square matrix T(z) whose 
elements are given by ©. An important feature of 
this matrix is that the column vector v with elements 
Vki = Ski is always left eigenvector of F(z ) with eigen- 



value z 1 , i.e., J2k T 



This equal- 



ity ensures the conservation of the trace of the density 
matrix p, and follows simply from ^2 k (k\Tr£(. ■ -)\k) = 
Tr(. . .). The matrix T(z) can thus be written as T(z) = 
z _1 u(z)v T + T'(z) where v T r'(z) = and v T u(z) = 1. 
The column vector u(z) is right eigenvector of T(z) with 
eigenvalue z~ l . Provided it has no pole on the real axis, 
the corresponding term of T(z) can be analytically con- 
tinued in the lower half plane and gives a constant con- 
tribution to the time-evolved density matrix Since 
v T r = X)fe r fefe = 1 f° r an y density matrix p(0), this 
contribution does not depend on the initial state of S. 
In summary, the steady state of the open system S is 
i Uki\k)(l\ where u k i are the elements of the column 
vector u determined by 



lim {ir/T(ir))} u = u. 

»7->0+ 



(7) 



Note that the condition [p£,Hs] = was not used to 
derive this equation. 



2. Weak coupling limit 

To determine the steady state of S in the limit of weak 
coupling to £ , we first expand the matrix elements (JB]) in 
powers of the Liouvillian £,i nt . . . = [Hi ntl . . .]. We obtain 

iW*) = — — + * w(z) } (8) 



z - e k + ei 



z - e k > + ei' 



up to second order, where Jkl,k'l'(z) can be ex- 
pressed in terms of the correlation functions C$,, (t) = 
Tr[ps exp(itHs)hki exp(— itHsjhk'v] of the environment 
operators h kl = (k\H int \l) = hj k , as 

lu,m>{z) = dte^{e u —Cl l k ,(t) + e^'CHi-t) 

£ [5ive^C%,{t) + Suy^C^i-t)] }. (9) 
j 



In this expression, we have used the notation uj k i — 
t k — t\. The stationarity of the correlation functions C¥,,, 
stems directly from the steady environment assumption 
[P£, Hs] — 0. For the Hamiltonian ([1]) and with the defi- 
nitions ([2]), /in = -ha = vi + v 2 , h 22 = -h 33 = vi —V2, 
hi3 = h 2A = wi, hi2 = h 3i = w 2 and h u = h 23 = 0. 

In the absence of interaction between S and £, the 
eigenvalue problem ([7]) reduces to {tk—£i)uki = 0. Conse- 
quently, the only matrix elements u k i with nonvanishing 
zeroth-order approximations are that for which e k = 
Thus, if the energy spectrum {e k } is nondegenerate, the 
corresponding steady density matrix is diagonal in the 
basis {\k}}. In the opposite case, there can exist co- 
herences between states \k) of equal energy. The matrix 
elements u k i to zeroth order, are determined by the equa- 
tions 



E 

c t /=e 



lkl,k'l'(iO + )Uk'l' = 



(10) 



where k and I satisfy e k = The remaining coherences 
u k \ are at least of first order in H int . By writing explic- 
itly the coefficients jki.k'k' (*0 + ), it can be shown that, 
for an environment in thermal equilibrium, i.e., pg cx 
cxp(— Hg/T), the thermal state u k i oc 5 k i exp(— efc/T) is 
solution of (|10|) . even in the presence of degeneracy in 
the spectrum of H$, see Appendix. 



B. Different splitting energies 

For unequal nonzero Ai and A2, the spectrum of the 
Hamiltonian Hg = — is non degenerate. 

The zeroth-order steady state of the TLS is thus a statis- 
tical mixture of the states $2$ . The relation (fTU)) becomes 





li_ 

_ 7 2 _ 7i 



7i -7 2 7i + +lt\ 
7i + + % -72 
■ 7 2 + 





7 2 + 



-7i + 



\ 


(Pl\ 




P2 




P3 


) 





= 



(11) 



where p k = u kk . The elements of the above matrix can 
be written as 



= 2t:Y,Pa/b\(B\w 1 \A)\ 2 S(E a - E B + A,) (12) 



AM 



where Ea and \A) denote the eigenenergies and eigen- 
states of Hg, and Pa are the eigenvalues of pg. The 
coefficients and 7" are the Fermi golden rule rates of 
the TLS i [23. 

The solution of (fTTj) leads to a product steady state 
p = pi <8> p 2 where 

Pi = (t+%r 1 [7t\+)iii+\+%\-)ii{-\]' as) 

The two TLS are uncorrelated, to lowest order in H int , 
when their splitting energies are different from each 
other. Moreover, the steady state pi of TLS i is the 
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same in the presence or absence of the other TLS. In 
the special case A 2 = 0, the zeroth-order coherences U12 
and it 34 are a priori different from zero since £1 = 62 and 
£3 = £4. But, for an environment £ consisting of several 
heat baths, it is shown in the Appendix that p = pi® 7/4 
where pi is given by ([T3|) and I is the 2x2 identity ma- 
trix, is steady state. 

C. Identical splitting energies 

For Ai = A 2 = A ^ 0, the states |2) and |3) have 
the same energy £2 = £3 = 0. The other energies are 
£4 = — £1 = A. Here, equation (fTQ|) takes the form 

(-3- <"> 

where c = M23 = W32, 7 is the 4x4 matrix given in (jlip 
and p T = (pi P2 P3 Pa). The coefficient a reads as 

a = \ (tT + if + 12 + 7 2 + )+E Pa\(A\v\B)\ 2 5(lo ab ) 

1 A,B 

+ [l(^K|5)| 2 - \(A\w 2 \B)\ 2 ] Pa + Pb (15) 

where w = 2y / 7r(vi — U2) and = Ea — -Eg- The 

elements of j3 T — {(3\ /?2 /?3 AO are given by 

1/4 = 27r5^P A/B <A|t« 1 |B)<fl|«;J|A)*(w i 4B-A) ) (16) 

A,B 

P2 = -(J3i+Pi)/2+ij3 and/3 3 = -(fit + J 8 4 )/2-ij8, where 
/9 = EAB(^l w 2|5)(5kil^X^-Ps)/(^i3 + A). For 
two identical two-level atoms coupled to the electomag- 
netic vacuum, /34 and /3 are, respectively, the collective 
decay rate and the dipole-dipole interaction energy of the 
atoms [Hill. 

It is instructive, for the following, to relate the co- 
efficients (fTBf to Fermi golden rule rates. Instead of 
analysing the influence of £ on the two TLS in the basis 
of product states ©, the basis made up of the states |1), 
1 4) and the entangled Bell states 

= = -^(l+>i|->2 ± |-}i|+> 2 ), (17) 

can be used. Both bases correspond to the same en- 
ergy spectrum {±A, 0}. The Fermi golden rule rates 
for the downward transitions |4) — > l?/^) are given by 
27Tj2 AB P A \(B\(^\H mt \4)\A)\ 2 S(u JBA - A) = (7+ + 

72 )/2 ± Re/^. This last expression is also valid for 
\tp ) — » |1). For the upward transitions |1) — > and 
1^) |4), the rates are (7^ + 72~)/2 ± Re/3i. 

Equation (TT4")) can be solved by diagonalizing 7. The 
eigenvalues of 7 are Aq = 0, Ai = 7T + 7i j ^2 = 7T + 



J2 and A4 = Ai + A2. We denote by ip n and 4> n the 
corresponding right and left eigenvectors. Since x]jq is the 
only right eigenvector for which the sum of its elements 
does not vanish, p = tp + J2 n >o ^^n^ni ! 3 * + c*/ 3 ); 
and the coherence c is solution of 

ac - (3 T *p - ]T A-^V^O^c + /3c*) = 0. (18) 

n>0 

We will see in the next section that c can be nonzero 
and lead to stationary entanglement of the TLS. In the 
special case Ai = A2 = 0, it can be shown that c = 
and pk = 1/4 are solutions of (fT"4"|) , if £ is made up of 
heat baths, see Appendix. 

The treatment of section ILU Bl applies when the differ- 
ence 6 = Ai — A2 is large enough that it can be considered 
finite in the expansion in terms of the interaction Hamil- 
tonian H int . In (fl"4"|) . this difference is exactly zero. A 
possible approach to understand the influence of a small 
6, consists in expanding the coefficients © both in Hi nt 
and S. This gives equation (|14[) with a + id in place of a, 
which reduces to (Q~4|) for 6 much smaller than the other 
matrix elements, and leads to the uncorrelated state (fT5)l 
with Ai = A 2 , in the opposite limit. 

IV. MULTIPLE HEAT BATHS ENVIRONMENT 

In this section, we consider an environment £ made up 
of several heat baths, as sketched in FigjTJ each consist- 
ing of an infinite number of harmonic degrees of freedom 
which are coupled linearly to the TLS. In other words, the 
spin-boson model [23[ , which appropriately describes var- 
ious physical environments fl9l | , is generalized to two TLS 
and several heat reservoirs. We show that two bosonic 
baths can induce stationary entanglement of two identi- 
cal TLS. 



A. Environment model 

We write the Hamiltonian of £ as He — J2 n ^Sn where 
n runs over the heat baths and 

Hen = y"^ W ga»o a "<z- ( 19 ) 

8 

In this expression, the sum runs over the harmonic modes 
of the bath n. The annihilation operators a nq satisfy the 
bosonic commutation relation [a n q,a n , ql ] — S nn '6 qq '. For 
the coupling operators, we consider 

and a similar expression for i?j . The coupling parameters 
knq are assumed to be real. The environment is initially 
in the state pe oc rj exp(— He n /T n ) where T n is the 
temperature of bath n. 
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Here, the rates (fT2"|) can be written as 

^^EldtUr (21) 

n 

where J,^' = 2nJ2 q [knq} 2 S(w nq — Aj), and the coeffi- 
cients (|16|) . which are relevant only in the case Aj = A, 
are given by similar expressions with J„ replaced by 
K n = 2irJ2 q knqknqS(uj nq - A). Clearly, J$ is neces- 
sarily positive but not K n , and \K n \ 2 < J n ^ J n ' ■ The 

(i) 

main difference between J rl and -ftT„ is that the former 
depends only on the coupling of TLS i to £, whereas the 
latter is determined by both coupling operators w\ and 
u>2- An important physical parameter that controls the 
ratio \K n \ 2 /JnJn is the distance d n between the TLS 
coupling points to bath n. This ratio reaches its maxi- 
mum value of 1 when the two TLS interact in exactly the 
same way with bath n, which necessarily means d n = 
13]. This spatial dependence is discussed more fully at 
the end of section IfV C 21 

Finally, we comment on the second term in (|15[) , which 
plays a role in the following. It can be cast into the 
form ^2 n J Q dui5{ui)L n {ui) / tanh(w/2T„) where the spec- 

(i) 

tral functions L n are denned similarly to Jn with u 
in place of Aj. The function L n vanishes for frequen- 
cies oj higher than a cut-off frequency 19]. Its low- 
frequency behavior leads to various possibilities. First, 
a is finite only if, for any n, L n /uj does not diverge for 
uj — » 0. If, in this limit, this ratio goes to zero for any 
n, then the second term in (|15p vanishes. This term 
reads as a n T n + a n 'T n i + . . . for Ohmic spectral densities 
L m ~ uj, m = n, n' , ... 19]. For a bath consisting of a 
Z?-dimensional continuous field, L n ~ oj d for uj — > 0, and 
hence is Ohmic for D = 1. However, note that, whereas 
Jn and K n are determined by the transverse coupling 
operators Wi, the functions L n depend on the longitudi- 
nal coupling. Consequently, L n can in principle be made 
as small as we wish, irrespective of the transverse cou- 
pling strength. 



B. Steady state for identical two-level systems 

From now on, we consider the case of identical TLS 
splitting energies Aj = A, for which, as seen above, sta- 
tionary TLS entanglement may exist. We further assume 
that the two TLS are coupled identically to the heat 
baths, i.e., Jn^ = J n = J n - This can hold for w\ ^ w% 
if the two TLS are connected to different points of bath 
n. As a consequence of these assumptions, 7^ = \ , see 
(|21[) . To simplify the following expressions, we introduce 
the notations : 

li = it = 1 , 7i~ = = IV (22) 
f3 4 = (3 , ft = /V , a = 7 (1 + V + 



where j3, rf and £ are real for the coupling operators ([20]) 
and with the above assumptions. As discussed above, £ 
is determined by the longitudinal coupling, whereas all 
the other parameters are related to the lateral coupling 
operators Wi. The coefficient $, denned right after (|16p . 
is also real and does not contribute to the TLS steady 
state. 

Under the assumption of real /3, r[ , £ and /?, we find, 
from (fl"8|) , a real coherence 

c=j(v'-v)(i-v 2 ){^i-v) 2 + f(i+v) 3 (^+v+0 

- (l + 7/)(l + 377 2 + 37/ + 7?y)} . (23) 

The populations of the TLS steady state p can be written 
in terms of c as 

( p 7/ 3 )^(-9 + ^"' + "' ) (^ 1 ) 

+ (7/77 + 7/ - V - 1) ( l-7?j I (24) 

where z — 1 + r\. Note that c = and p is uncorrelated 
for rf = r) or 77 = 1. When this last equality is satisfied, p 
is proportional to the identity matrix, as expected from 
the case A = 0, see Appendix. The equality r/' = r/ 
holds, for instance, when £ is in thermal equilibrium. 
The denominator in (|23|) vanishes for 7 = /?, 77' = 77 and 
£ = 0. These three conditions are fulfilled for W2 = w% 
and vi = v\. There is not a unique steady state when 
the two TLS interact with £ in exactly the same way [l3| . 
We also remark that, since c is real and p 2 = p 3 , the TLS 
steady state can be written as p = p 1 |l)(l| +p 4 |4)(4| + 
(P2 +c)|-0 + )(V> + | + ipi ~ c)\ip~){ip~\ with the Bell states 
lip^) given by pT)l . We wil see below that, though 
and \tp~) have the same energy £2 = £3 = 0, there exists 
a parameter regime in which 754 = and pi = |c|, and p 
is hence entangled. 

C. Entanglement induced by two heat baths 

We now study the entanglement of p for an environ- 
ment £ that consists of two heat baths of temperatures 
T\ and T 2 . The steady state p is entangled if and only if 
its partial transpose p T = Y, Pk\k)(k\+c(\l)(4\ + |4)(1|) 
has negative eigenvalues [24l 1251] . The eigenvalues of p T 
arep 2 = p 3 and A± = (p 1 +p A )/2±[( Pl -p 4 ) 2 +4c 2 }^ 2 /2. 
Clearly, only A_ can be negative. 

1. Low-temperature entanglement region 

As an interesting example, we consider the case £ = 
and K\ = 0. This last condition means that the indi- 
rect interaction between the TLS is mediated only by 
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T| (units of A) K-i/Ji 



FIG. 2: Entanglement region in the (Ti, T2) plane for K\ — 0. 
The TLS steady state is entangled for temperatures below the 
drawn line. The solid lines correspond to K2 — J2, £ = 0, 
and J2/J1 ~ 5, 10 and 50. The size of the entanglement 
region increases with J2. For the dashed and dotted lines, the 
coupling parameters are f = 0, and, respectively, K2 = 0.72 J2 
and J 2 = 500 Ji, and K 2 = O.8J2 and J 2 = 150 Ji. The short- 
dashed line is obtained for an Ohmic £ = O.OlTi/A and for 
K 2 = J2 =50Ji. 



bath 2. With this value of K\, the results discussed here 
hold also for the three bath setup depicted in FigfT]when 
T3 = T\. We find that there can be a low-temperature 
region, determined by J2/J1 and if 2/ Ji, in which p is 
entangled, see FigJ^l We remark that the line delimiting 
this entanglement region in the (Ti,T 2 ) plane, is tangent 
to the equilibrium line T± = T 2 for Tj.,22 — > 0, and is 
essentially vertical at its other end for T2 -C A. These 
two behaviors come from the fact that the temperatures 
contribute to p only via Boltzmann factors exp(— A/T„). 



Analytical results can be obtained by expanding the 
eigenvalue A_ to lowest order in these factors. It as- 
sumes negative values in the vicinity of Ji = T% = 0, for 
\K 2 \ > J2/V2 and J 2 > Ji[V2\K 2 \/J 2 - These re- 

quirements are the same in the Ohmic case discussed at 
end of IIV A[ for which £ = £iTi+!; 2 T 2 vanishes in the lim- 
its Ti,T 2 — > 0. For given coupling parameters satisfying 
the above conditions, p is not entangled if the tempera- 
tures Ti and Xb are too high. However, for £ = 0, the 
maximum possible value of T\ is proportional to AJ 2 / J\ 
in the large J 2 limit, see Figj2j Consequently, in this 
case, entangled states exist for any temperature T\. For 
T2, in contrast, our numerical results suggest that the 
steady state is always separable for T 2 greater than a 
value of about 0.567A. Entangled states can be observed 
close to this temperature in the limit of large J 2 . For 
£ = £iTi + £2X2, p is necessarily separable for T\ higher 
than a temperature that diverges for £1 — > 0. 



FIG. 3: Region of the parameter plane [K\j J\,Kij J-z), K2 > 
0, where entangled steady states can be found, for J2/J1 = 
2.4, 2.7, 3, 5, 10, 100 and 1000. The entanglement region is 
above the drawn line. Its size increases with J2, to a maximum 
asymptotic value which is practically reached for J2 = 1000 Ji . 



2. Requirements on the characteristics of the environment 

Stationary entanglement can also be obtained for K\ ^ 
0. Since p is obviously invariant under the bath per- 
mutation (Ji,Ki) <H> {J 2 ,K 2 ), it is enough to consider 
J 2 > J\. In this case, it can be shown that there exist 
entangled steady states in the vicinity of T\ — T 2 = if 



(Ji + -h) - (Ki + K 2 f - \Kj. + K 2 



J-2 



K 2 - Ki-^- < 0. 

(25) 

This condition remains the same if the signs of both Ki 
and K 2 are changed. Figure [3] shows, for K 2 > 0, the 
coupling parameter region where p can be found entan- 
gled. The following interesting conclusions can be drawn 
from these results. There is a particular value of J 2 /J\ 
below which p is separable. In other words, the couplings 
to the two heat baths must differ enough from each other 
in order to observe stationary entanglement. For given 
J 2 / J 1 and K\l J\ such that entangled steady states exist, 
these states are obtained for \K 2 \/J 2 not too far from 1. 

As mentioned above, the ratio \K 2 \/J 2 depends essen- 
tially on the distance d 2 between the two points of bath 2 
where the TLS are connected. More precisely, it is deter- 
mined by a dimensionless parameter d = Ad 2 /v where 
v is a characteristic field velocity of bath 2. The ratio 
I K 2 \/J 2 is small for large d. This imposes limitations 
on A and on the temperature X2 to obtain an entangled 
steady state. A distance d 2 of l//m and a low field ve- 
locity v of 10 3 m.s -1 , which is the order of magnitude of 
the sound velocity in solids, give a temperature of about 
10 mK, which is an experimentally accessible value. An- 
other important characteristic of bath 2 is its dimension- 
ality D. For example, for a continuous free field, K 2 /J 2 
is equal to sm(d)/d for D = 3, Jo(d) where Jo is the ze- 
roth order Bessel function of the first kind, for D = 2, 
and cos(g?) for D = 1. Thus, in this last case, stationary 
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0.08 0.16 

Ti (units ofAJi /J 2 ) 



FIG. 4: Negativity as a function of Ti in units of AJ1/J2 for 
K 2 = J 2 = 100 Ji (short-dashed line), K 2 = J 2 = 1000 Ji 
(full line), #2 = 0.95 J 2 = 9.5 Ji (dash-dotted line) and K 2 = 
0.95J2 = 95 Ji (dashed line). The dotted lines correspond to 
the large J 2 approximation discussed in the text. The inset 
shows the populations of the ground and singlet states as 
functions of Ti in units of AJi/J 2 for K 2 = 0.95J 2 = 95 Ji, 
K 2 = J 2 = 100 Ji and K 2 = J 2 = 1000 Ji. The two other 
populations are small. The other parameters are £ = 0, iii = 

and T 2 = 0.1A. 

entanglement can be obtained for large distances d 2 and 
the limitations discussed above do not apply. 

3. Maximum attainable entanglement 

Finally, we present quantitative results for the entan- 
glement of the steady state p. As a measure of entangle- 
ment, we use the negativity M(o) = (||p r ||i — l)/2 where 
||.||i denotes the trace norm [!, [26j]. Negativity ranges 
from for separable states to 1/2 for maximally entan- 
gled states. Here, it is equal to — A_ when this eigenvalue 
is negative, and to otherwise. The maximum value of 
N that we have found, is reached for the coupling param- 
eters K\ = 0, £ = and J 2 3> Ji, and the temperatures 
T 2 < A and T x > A, see Figg] In this regime, the TLS 
steady state is given by 

pi = a [(1 + 6>) 2 (1 + 26) - k 2 ] , p 2 = 0-6/(1 + + 26») 
Pi = a9 2 (l + 29) , c=-okQ (26) 

where a = [(1 + 26>) 3 - k 2 ]" 1 , k = K 2 /J 2 and 6» = 
(Ti/ A)(Ji/ J 2 ). For |k| 7^ 1, as 9 increases from zero to 
infinity, the ground state population p\ decreases from 

1 to 1/4, p2 = P3 and P4 increases from to 1/4, and 
|c| increases from zero to a maximum and then decays 
back to zero. The low 9 behavior is very different for 
\k\ = 1. In this case, c is finite in the limit 9 —> 0. 
For K 2 — ±J 2 and a temperature A <C T\ <C AJ 2 /Ji, 
we find p = (2/3)|+)n(+| <g> |+) 22 (+| + (1/3)|^)(^,T| 
with the Bell states |-0 =F ) given by flT"7l) . and a negativity 
A/" = (-\/5 — 2)/6 ~ 0.04. The same entangled state can be 



reached for K\ 7^ 0, as it will be clear from the discussion 
below. Our numerical results suggest that finite values of 
A" correspond generally to states p such that essentially 
only the ground state and one of the Bell states (fT7|) are 
populated, see inset of Figj4] 

To better understand the above results, it is interesting 
to consider the rates discussed after (fT?]) . For T 2 <C A 
and T\ 3> A, the rates of the upward transitions |1) — > 
IV^) and |^ ± ) -> 1 4) are r± p = ( J x ± K^Tx/A, and that 
of the downward transitions |4) — > |i/' ± ) and j^^) — > |1) 
are J 2 ± K 2 + r± p . For X 2 = J 2 and T x < AJ2/J1, 
the rate of |4) — > and |^ _ ) — > |1) is equal to r~ p , 

whereas that of |4) -> \ip + ) and -> |1) is 2J 2 > r^ p . 
Consequently, the states |4) and \4> + ) are essentially not 
populated, and the transition rate from the state IVO 
to the ground state |1) is effectively twice that of the 
reverse transition, leading to a factor of two between the 
two corresponding populations. The situation is similar 
fo K 2 = — J 2 . If \K 2 \ is too far from J 2 or if T\ is too 
high, the values of the different rates are comparable and 
so are the populations of the states |1), |i/' ± ) and |4), and 
hence p is separable. 

V. CONCLUSION 

In summary, we have studied a system of two indepen- 
dent TLS weakly coupled to a stationary non-equilibrium 
environment. Considering first a general open system, 
we have determined their steady state. Without speci- 
fiying any further the surroundings of the TLS, it can 
be shown that their steady state is uncorrelated if their 
splitting energies are different from each other. More- 
over, the state of each TLS is the same wether or not the 
other TLS is present. Consequently, in this case, a finite 
strength of the coupling to the environment is required to 
possibly generate stationary TLS entanglement. In the 
opposite case of identical splitting energies, on the con- 
trary, stationary correlations between the TLS can exist 
for extremely weak coupling to the environment. 

To determine wether these correlations can be quan- 
tum, we have considered the case of an environment con- 
sisting of several bosonic heat baths at different temper- 
atures. We have shown that, for TLS coupled similarly 
to two baths, there are temperatures and coupling pa- 
rameters for which the TLS steady state is entangled. 
An important requirement is that, for at least one bath, 
the points to which the TLS are connected must be close 
enough to each other. However, this condition can be 
relaxed when one of the bath is one-dimensional. In this 
case, the TLS can be as far apart as we like. There are 
also requirements on the baths temperatures. Essentially, 
one of them must be sufficiently low, of the order of the 
TLS splitting energy. Depending on the characteristics 
of the coupling, the other temperature can be unlimited. 

We have found a parameter regime where the TLS 
steady state is a statistical mixture of the product ground 
state and of the entangled singlet state with weights 2/3 
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and 1/3, respectively. This mixed state is entangled and 
the corresponding negativity is about 0.04 which is the 
largest value we have obtained. Interestingly, this regime 
can be fully understood in terms of Fermi golden rule 
transitions between appropriate states. To conclude, our 
results show that a relatively simple non-equilibrium en- 
vironment can lead to stationary entanglement of two 
TLS, but certainly do not exhaust all the possible effects 
of stationary non-equilibrium surroundings on quantum 
correlations. Larger entanglement of independent TLS, 
as measured by negativity for instance, may be achievable 
with other environments or for TLS coupled differently 
to the environment. Further studies in these directions 
would be of interest. 



k = I. For k I, we start by showing that j 2 = 7^ 
which implies pi = p 2 and = p±, see l[i~3"|) . The differ- 
ence of these rates reads as 

7 2 + - 2itJ2(Pa -Pb)\{B\w 2 \A)\ 2 5(uj ab ). (A2) 

A,B 

For the kind of environment considered, Hg = Hg n 
and hence its eigenstates and eigenenergies can be writ- 
ten as \A) = U n \ A(n) ) and E A = E n E AM- The pop- 
ulations Pa factorise as Pa oc J7„ exp(— E A { n )/T n ). The 
TLS are coupled to each bath thus w% — w 2n ■ Con- 
sequently, the difference (|A2[) satisfies 



Appendix A: Special uncorrelated states 

In this appendix, our purpose is to show that, for some 
special cases, the solution of (flU)) is of the form Uki — 
PkSki- This is the case if the sums 

J2"fki.k>k'(i0 + )Pk' =tt PA{A\h kk ,\B){B\h kn \A) 

k' A,B.k' 

^ ( 1 Pb \ xi \ 1 * VI- Vk 

X [Pk +Pi - 2-^rPk> o(uiab - Wfc/fcJ H 

\ Pa J 7T uj A b - uj k 'k 

(Al) 

vanish for k and I such that tk — 

For an environment in thermal equilibrium, i.e., 
Pa oc exp(— Ea/T) where T is its temperature, Uki oc 
exp(—ek/T)6ki satisfies (fiT)]) since, in the sums (|A1[) . 
p k =pi and P B Pk'/PAPk = exp[-(oj B A +u k 'k)/T}. This 
proof applies to any system S. 

We now consider the case of zero splitting energy A2 
and of an environment £ that consists of heat baths at 
different temperatures T n . First, the populations pk ob- 
tained in section IIIIBI ensure the vanishing of (|A1|) for 



x (e- E A^I^ - e^wA) S(E A(n) - E BM ) (A3) 

and hence vanishes. For A2 = 0, the sum (IA1|) must be 
zero for (k,l) = (1,2), (2,1), (3,4) and (4,3). For these 
cases, the equalities p\ = p 2 and p^ = p±, shown above, 
lead to 

J2w(W + )Pk> =2tt (A\h kk ,\B){B\h kn \A) 

k' A,B,k> 

X (p k PA - P B Pk')S{uJAB - Wfc/fe). (A4) 

For the Hamiltonian |T]), hn = h 2 3 = and hence the 
only terms that contribute to the above sum are such 
that p^ = Pk and u>kk' = 0. Thus, it vanishes for the 
same reasons as (|A2[) does. 

For Ai = A2 = 0, the sum (|A1[) must be zero 
for any (k,l). In this case, for an environment that 
consists of heat baths, Uki — S k i/4 satisfies (fTOjl since 
E k >lki,k>k'(i0 + ) = 0. 
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